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Abstract 

We argue that claims about magnetic field dependence of the magnetic field pene- 
tration depth A, which were made on the basis of muon-spin-rotation (^uSR) studies 
of some superconductors, originate from insufficient accuracy of theoretical models 
employed for the data analysis. We also reanalyze some of already published exper- 
imental data and demonstrate that numerical calculations of Brandt [E.H. Brandt, 
Phys. Rev. B 68, 54506 (2003)] may serve as a reliable and powerful tool for the 
analysis of /LiSR data collected in experiments with conventional superconductors. 
Furthermore, one can use this approach in order to distinguish between conventional 
and unconventional superconductors. It is unfortunate that these calculations have 
practically never been employed for the analysis of ^SR data. 
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1 Introduction 



Muon-spin-rotation (/iSR) experiments in the mixed state of type-II supercon- 
ductors provide unique information about superconducting properties of the 
investigated sample. An important advantage of this method is that muons 
probe the bulk of the sample and therefore, the results are not distorted by 
possible imperfections of the sample surface. At the same time, in order to 
extract quantitative results from /iSR measurements, a detailed model of the 
magnetic field distribution in the mixed state is needed. As well as we are 
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aware, only the Ginzburg-Landau (GL) theory [T] of the Abrikosov vortex lat- 
tice [2j is developed to such a level [31IHI5] . As was recently demonstrated, if an 
adequate model is available, not only the magnetic field penetration depth A 
but also the upper critical field Hc2 can be found from /iSR data collected in 
different applied magnetic fields [6]. It has to be remembered, however, that 
theoretical calculations of Refs. [3|Hl5] are related to superconductors with one 
and isotropic energy gap only. This is why, this kind of analysis should be used 
with extreme caution in the case of unconventional superconductors, in which 
the applicability of theoretical models is not obvious. 



We also point out a very interesting and promising approach which was de- 
veloped in Refs. [TllSjl^piI] . In these works a microscopic theory was used for 
calculation of the mixed state parameters. An important advantage of this ap- 
proach is that the results are not limited to conventional superconductors and 
it can be used at temperatures well below T^. both for s- and (i-wave pairing. 



In recent years, /iSR measurements were widely used for studying of different 
unconventional superconductors such as high-T^ materials, MgB2 and oth- 
ers. Some very interesting results were obtained. It was demonstrated that in 
some cases the magnetic field penetration depth A and the superconducting co- 
herence length ^, evaluated from /^SR measurements, depend on the applied 
magnetic field (see, e.g., p!T]|12lll3|14|15|16|17lll8lll9] ). This result, however, 
contradicts the GL theory, which was used as a basis for the data analysis. 
This contradiction is a clear sign that the corresponding models are not ade- 
quate for describing the magnetic field distribution in the mixed state of these 
compounds and rises the question about physical meanings of \{H) and i{H) 
obtained in such a way. As we argue below, magnetic field dependences of A 
and ^ cannot be obtained from /iSR experiments if the conventional GL the- 
ory or the London model were employed for the analysis of experimental data. 
Moreover, because in the mixed state the superconducting order parameter is 
not spatially uniform, there is no reasonable way to define either A or ^. In 
other words, the physical meanings of magnetic field dependences of A and 
^, evaluated from /iSR data, are quite different from traditional definitions of 
these two lengths. This circumstance was recognized in Refs. [2(111211122] where 
it was pointed out that A(iJ), evaluated in such a way, represents some fit- 
parameter rather than the magnetic field penetration depth. We underline 
that the same should also be addressed to ^ {H) dependences. In the following 
section, in order to avoid confusion, we shall use Aq and to denote values A 
and ^ for H ^ 0. 
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2 Conventional superconductors 



Superconductors with s-pairing and one energy gap we sliall consider as con- 
ventional, independent of their pairing mechanism. Because the GL theory is 
traditionally used for analyses of /iSR data, we limit our consideration to this 
theory. 

The magnetic field penetration depth Aq together with the zero-field coher- 
ence length ^0 represent two fundamental lengths of the GL theory. If their 
values for some particular temperature T are known, one can calculate the 
GL parameter 

k{T) = Xo{T)/UT), (1) 



the thermodynamic critical magnetic field 

HJT) = , (2) 



the upper critical field 

H,,{T) = V2kH,{T) = -J^, (3) 



the lower critical field 



with a{K) = 0.49693 + exp [-0.41477 - 0.775 In k - 0.1303(lnfi:)2] [5] Further- 
more, in the case of conventional superconductors, any characteristics of the 
sample for any value of an applied magnetic field may also be calculated and 
expressed via Aq and ^o- Very detailed numerical calculations of different pa- 
rameters of the mixed state for a very wide range of k and for magnetic fields 
ranging from Hd to Hc2 are presented in Ref. [5]. 

Muons probe the distribution of the magnetic induction in the sample. In high- 
K superconductors and low magnetic inductions B, contributions of vortex 
cores can be neglected (London limit) and the distribution of the magnetic 
induction around a single vortex line may be written as 

Bir) = A^o(r/Ao), (5) 
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where r is the distance from the vortex center, $o is the magnetic flux quantum 
and Kq is the modifled Bessel function. Because Eq. (5) is obtained from the 
London theory, it gives an unphysical divergence of S at r = 0. In order 
to improve Eq. (5), an appropriate cutoff has to be introduced j231l24f25] . It 
should be remembered, however, that the results of Refs. [23|24|25j can be 
considered as sufficiently accurate in low magnetic fields H <^ Hc2 only. If 
this condition is not satisfied, numerical solution of the GL equations must 
be used for a reliable analysis /iSR data. The magnetic induction distribution 
may be calculated as a linear superposition of inductions created by different 
vortices (see, for instance, Ref. [23]). 

By measuring muon relaxation rates, one obtains the distribution of the mag- 
netic induction P{B) experimentally, which allows to calculate the variance 
of the magnetic induction 

= {B^{r)-W) , (6) 

where = J ■ ■ ■ d^r means spatial averaging over superconductor of 

volume V . If the distribution of the magnetic induction around vortices is 
known, a can also be calculated theoretically. According to [5] 

a = F{K,B/B,2)/Xl (7) 

where the parameter F depends on k and B/ Bc2- If the value of F is known, 
Ao may straightforwardly be evaluated. In the case of k 3> 1 and 6 <^ 1, 
F ^ 0.061$o- In other situations, reliable results can be obtained from Ref. 
[5]. Eq. (7) may also be written as o" = {27rHc2/^o)F{K,B/Bc2)/K'^. This 
representation may be convenient if evaluation of k is preferable. 

While the zero-field value of A enters the theory, the actual magnetic field 
penetration depth is field dependent. According to the original Ginzburg and 
Landau publication [1], if the magnetic field is parallel to the sample surface, 

X{H) = Ao [1 + f{'^)H/H,] . (8) 

The function /(k) grows monotonically with k in such a way that for k <^ 1 
/(k) ~ k/4:V2 and f{oo) = 0.125 [1]. Taking into account Eq. (1), we see 
that even the magnetic field dependence of A may be expressed via Aq and 
^0- The \{H) dependence arises due to suppression of the order parameter 
\ip\ by the applied magnetic field. In bulk type-II superconductors, Eq. (8) 
is applicable in the Meissner state only, i.e., in magnetic fields H < Hd- If 
H > Hci, the magnetic field penetrates into the bulk of the sample forming a 
lattice of Abrikosov vortices. 
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Fig. 1. F multipied by (1 — 0.069/k^) as a function of B/Bc2 according to [5j. The 
dashed hne shows the F{B/Bc2) according to interpolation formula proposed in [3j. 
The horizontal line corresponds to -F = 0.061<l>o- The inset shows the same curves 
on linear scales. 

If spatial variations of the order parameter can be neglected, the magnetic 
induction decays exponentially on the fiat surface of the sample. In the case 
of cylindrical geometry (around vortices), the same decay is described by the 
Bessel function (see Eq. (5)). Considering /xSR experiments in the mixed sate 
of type-II superconductors, we can use Eq. (5) if the total volume of vortex 
cores is negligibly small compared to the volume of the sample, i.e., k ^ 1 
and B <^ Bc2- If one or both of these conditions is not satisfied, spatial 
variations of the order parameter have to be taken into account. Numerical 
calculations of Brandt are shown in Fig. 1. If k 3> 1, the function F{B/Bc2) 
is practically independent of k. As may be seen in Fig. 1, F remains magnetic 
field dependent even for very small values of B/Bc2- The maximum on the 
F{B/Bc2) dependence at B/Bc2 ~ is the obvious consequence of the 

fact that a vanishes at B/Bc2 — >■ 0. At very low magnetic inductions F is 



It must be clearly understood that, although the applied magnetic field influ- 
ences properties of the sample both in the Meissner and the mixed states, the 
physics of this influence is completely different. In the Meissner state, super- 
currents are induced in the surface layer of the sample. The density of these 
currents is proportional to H and they depress the order parameter, which 
leads to an increase of A (see Eq. (8)). Because the reduction of the order 
parameter is small (|A?/'| <^ I'ipl), one may still introduce the magnetic 
field penetration depth in its traditional way. 

In the mixed state, the situation is completely different. Because there are 



proportional to J{B/Bc2) ^5]. 
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no currents, which are proportional to H, the absolute value of the applied 
magnetic field is irrelevant. Only the distance between vortices given by the 
magnetic induction B is important. At magnetic inductions B < 0.1Bc2, over- 
lapping of vortex cores may be neglected and vortex properties are indepen- 
dent of the applied magnetic field [5]. Only the vortex density is changed. At 
higher magnetic inductions, vortex cores overlap and not only the vortex den- 
sity, but also properties of individual vortices are magnetic field dependent. 
Because the local value of the magnetic field penetration depth is inversely pro- 
portional to the modulus of the order parameter {tpl, there is no much sense 
to introduce any unique value of A corresponding to each particular magnetic 
field. The correct approach is to calculate some measurable quantities theo- 
retically and compare them with experimental results. In this way, however, 
only the zero-field value of A can be evaluated. If the value of A resulting from 
the analysis of experimental data depends on the applied magnetic field, it 
means that the theory, which was employed for the analysis, does not describe 
the actual experimental situation and the approach to the analysis should be 
reconsidered. 

If the magnetic field dependence of F is not taken into account or accounted 
for incorrectly, the analysis of muon depolarization rates would result in some 
effective Xeff, which is magnetic field dependent. The knowlege of \eff{H), 
however, does not represent any particular interest. This is why it is important 
to use reliable models of the mixed sate in order to evaluate Aq. 

As is well known, the GL theory is formally applicable at temperatures close 
to Tc only. This is why quantitative applicability of theoretical calculations to 
the analysis of experimental data at temperatures well below is not obvious. 
However, as it was recently demonstrated, the magnetic field dependence of 
a at T can be very well fitted by calculations of Brandt with two fit- 
parameters Ao and Hc2 [6]. Moreover, the value of Hc2, evaluated in such a way, 
coincides with the result of magnetization measurements. We consider this as a 
proof that the theoretical dependence calculated in framework of the GL 
theory can indeed be used for quantitative analysis of isothermal experimental 
data even at temperatures T <^ Tc. Below we reconsider several experimental 
/iSR studies and demonstrate that their results may perfectly be described by 
the conventional GL theory although the magnetic field dependence of A was 
claimed in some of the original publications. 



3 Unconventional superconductors 

It must be remembered that calculations of Brandt [5j, which we have dis- 
cussed above, are valid for conventional superconductors only. There are no 
reasons to believe that the vortex core structure should be the same in two- 
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gap superconductors or in superconductors with nodes in the order parameter. 
Furthermore, one may assume that the influence of the vortex core region on 
the distribution of the magnetic induction should be even stronger than in 
the case of conventional superconductors. This is why, if calculations of ^ or 
any other calculations based on the conventional GL theory are used for the 
analysis of /iSR data collected in different magnetic fields, it would produce 
an unphysical A (if) dependence. Although this result does not mean any spe- 
cial behavior of the magnetic field penetration depth, it should be considered 
as interesting. Indeed, if the conventional GL theory cannot describe the re- 
sults of fiSK experiments and all other possibilities for this disagreement are 
excluded0 one may conclude that this superconductor is unconventional. 

Superconductors with d-pairing as well as two-gap superconductors are more 
complex than conventional ones. For instance, two lengths and Aq are insuf- 
ficient for their characterization and some additional information is needed. At 
present, there is no experimentally proven theory of the mixed state in uncon- 
ventional superconductors. In this sense, magnetic field dependences of muon 
relaxation rates cannot be interpreted quantitatively without some additional 
assumptions. At the same time, one can try to obtain F{B/Bc2) experimen- 
tally in order to compare results for different superconducting materials. Un- 
fortunately, concerning high-T^ superconductors, the Hc2{T) curves are not 
yet reliably established. 

Interesting theoretical approaches for interpretation of the /iSR experiments in 
the case of d-pairing was developed in [71l8llM0ll26ll271l28]l29^ In [26ll271l28ll29] 
was convincingly argued that because of the nodes of the order parameter, 
the electrodynamics of the mixed state becomes nonlocal. This nonlocality 
effectively increases the vortex core radius and changes the distribution of the 
magnetic induction around vortices (see Fig. 6 of Ref. [2B])- If this effect is 
not taken into account, the magnetic field penetration depth evaluated from 
fiSR. experiments will be overestimated and magnetic field dependent. The 
distortion of the results is very clearly demonstrated in Fig. 4 of Ref. |28j . 
In order to correct the results, the function Ae//(-B)/Ao was introduced [28] . 
Using this function, which is an analog of F[B/Bc2), one can evaluate the 
magnetic field penetration depth Aq. 

At the same time, in high k superconductors and at low magnetic inductions, 
the total volume of vortex cores is small and contribution of vortex cores 
cannot considerably change the muon signal. In this case, one may use F = 
0.061$o for evaluation of Aq also in unconventional superconductors. Because 
A ~ resulting error is not expected to be big. This means that, if an 

experimental cr{H) dependence is available, extrapolation of cr(if) (or l/^/o") 



For instance, the traditional analysis cannot be used in the case of poly cry stalHne 
samples of anisotropic superconductors. 
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to H = gives more reliable values of Aq. 



4 Analysis of experimental results 

In this section, in order to simplify notation, we shall use A and C, without 
indexes, having in mind the magnetic field penetration depth and the super- 
conducting coherence length as they are introduced in the GL theory. 

4.1 RbOs206, Cd2Re20r, PrOs^Sbu- 

Experimental data for a polycrystalline sample of RbOs20e are shown 

in Fig. 2(a). Because RbOs206 is an isotropic superconductor, using of such 
samples is justified. This sample was investigated in |30] and experimental 
data were analyzed by employing of an interpolation formula proposed in 
[3j. Because this formula deviates significantly from more accurate numerical 
calculations (see Fig. 1), we reanalyze these data using calculations of Ref. 
[5]. As may be seen in Fig. 2(a), experimental data-points for H > 2 kOe 
can very well be fitted by the theoretical a{B/Bc2) curve. This fit results in 
H^2 = (67 ± 10) kOe and A = (220 ± 5) nm (the value of A = 260 nm was 
obtained in the original publiction). Because the value of ifc2 is obtained by 
the extrapolation of the cr{H) curve to a = 0, the corresponding error margins 
are large. It is important to emphasize that Hc2, evaluated in such a way, is in 
agreement with the Hc2{T) curve presented in [30]. This agreement together 
with sufficiently high quality of fitting strongly supports our analysis. 

We have chosen a high field part of the experimental cr{H) curve for the 
analysis because in higher magnetic fields F{B/Bc2) is independent of k (see 
Fig. 1). In principle, analyzing the low field part of the curve, the value of 
K, may straightforwardly be evaluated. This, however, is not always feasible. 
As was already mentioned, the correct parameter is not H but the magnetic 
induction B. The value of B determines intervortex distances and all other 
characteristics of the mixed state. In magnetic fields H ^ Hd, the difference 
[H — B) <^ H and one can use H instead of B. In low fields, however, the 
equilibrium value of B is considerably smaller than H and the actual difference 
{H — B) depends on pinning and on the demagnetizing factor of the sample. 
Furthermore, in low fields, the magnetic induction is nonuniform throughout 
the sample if its shape is not ellipsoidal. In polycrystalline samples, the sit- 
uation complicates even further. Indeed, in such samples, some vortices may 
go along intergrain boundaries, which can significantly infiuence the magnetic 
induction distribution. This is the reason that we do not speculate on the 
low- field behavior of the cr(if) curve. 
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Fig. 2. cr{H) data for three different superconducting compounds. The sohd hnes 
represent the theoretical a{B/Bc2) curve fitted to data-points. Only the data marked 
by closed symbols were used for fitting. The resulting values of A and Hc2 are 
indicated in the figure, (a) RbOs206 sample studied in pO]. The vertical dashed 
line indicates the value of Hd- (b) Cd2Re207 sample studied in [31]. The dashed 
line represents a linear approximation to a high field part of the theoretical cr{H) 
curve, (c) PrOs4Sbi2 studied in ^32j . 



Similar results for a Cd2Re207 sample studied in pi] are shown in Fig. 2(b). 
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For the reasons explained above, we disregard the lowest field data-point. 
Again in this case, data can be very well fitted with the GL theory, providing 
Hc2 = (5.75 ± 1) kOe in agreement with the original data (see [M])- The value 
of A = (830 ± 40) nm is also close to the result A = 750 nm of Ref. [HT]. 
We also note that approximation of experimental (t{H) data-ponts by a linear 
dependence, as it was done in [31] and some other publications, is unjustified. 
As may be seen in Figs. 1 and 2(b), the theoretical (y{B) curves are not at all 
linear. 



Fig. 2(c) presents cr(iJ) data for a heavy-fermion superconductor PrOs4Sbi2 
[22] • We do not discuss here different features of this rather unusual supercon- 
ductor but limit ourselves to one simple question whether the cr(if) depen- 
dence for this compound can be described by the conventional GL theory. 



As may be seen in Fig. 2(c) (see also Fig. 4 of Ref. [32]), the values of a{lkOe) 
and cr(2kOe) practically coincide. It was assumed in Ref. [32] that a change 
of vortex lattice symmetry or some other important changes of the vortex 
structure, which occur in magnetic fields above H = 1 kOe, may be responsible 
for such a behavior. This explanation seems to be plausible and we, as a 
precaution, do not use the highest field data-point in the analysis. 



The solid line in Fig. 2(c) represents the best fit of the theoretical cr{H) curve 
to the data collected in magnetic fields 0.2 kOe < H < 1 kOe. Quite amazingly, 
the resulting value of Hc2 = 21 kOe practically coincides with Hc2 = 22.2 kOe 
obtained in Ref. |33j from resistivity measurements. The value of A = (318±4) 
nm is also close to the result A = 290 nm of the original publication. 



Table 1 

RbOssOe (1.6K) CdaReaOr (O.IK) PrOs4Sbi2 (O.IK) 
Hc2 (kOe) 67 ± 10 5.75 ± 1 21.4 ± 4 

K 34 ±1 35 ±3 25.5 ±2.5 

A (nm) 220 ± 5 830 ± 40 318 ± 4 

^ (nm) 7 ±0.4 24 ± 3.5 12.5 ± 1.5 



The main characteristics of the superconducting compounds, resulting from 
our analysis of the /xSR data published in Refs. [5II|31ll32j . are listed in Ta- 
ble 1. We emphasize that all parameters were evaluated by fitting of cr(if) 
data-points with the theoretical a{B / Bc2) dependence calculated in [5]. In 
all cases, the values of Hc2 practically coincide with results of independent 
measurements. 
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4-2 CeRuiUB^. 



Calculation of cr(T, H) considered above is not the only way of analysis of /iSR 
experiments in the mixed state of type-II superconductors. A different method 
was employed in Refs. [TT|12|13lll4|15|16lll7|21|22j . In this approach, the dis- 
tribution of local fields (the Fourier transform of the muon precession signal) 
P{B) was directly analyzed by comparing with corresponding theoretical cal- 
culations. In real experiments, however, the P{B) line is usually different from 
theoretical predictions. This difference is expected. Indeed, the calculations are 
made for a perfect sample and for a perfect vortex lattice. All imperfections, 
which cannot be avoided in experiments, influence the P{B) curves. This is 
why, in order to approximate experimental data with theoretical calculations, 
some gaussian smearing factor is introduced. In such a way, satisfactory agree- 
ment between the theory and experiments can be achieved. This is justified 
if it is a priori known that the spatial distribution of the magnetic induction 
around vortex lines is in agreement with the theory, which is used for the cal- 
culations. If it is not the case, introducing of additional Gaussian relaxation 
may mask the disagreement and provide misleading results. 

For some type-II superconductors, the P{B) curves for different values of the 
applied magnetic field are available in the literature. This allows to calculate 
(j{H) and to employ the same kind of the analysis as was used above. Below we 
present the results of such analysis for single crystals of CeRu2 and vanadium. 

Fig. 3(a) shows a as a function of H for a CeRu2 sample experimentally 
investigated in Ref . pj] . The difference to the results displayed in Fig. 2 is that 
broadening of the P{B) line resulting from other sources of field inhomogeneity 
was not accounted for. In this case a, evaluated from /xSR experiments may 
be written as 

^ = + <, (9) 

where asc and abg are the mixed state and background contributions, respec- 
tively. As may be seen in [IB], (Jbg is not small and cannot be evaluated from 
the data presented in the publication with sufficient accuracy. This is the 
reason that we introduce a^g as an additional adjustable parameter. Because 
experimental data are insufficient for evaluation of A, iPa and ai,g together, 
we take the value of ifc2 from the original publicationo We also note that 
if = 40 kOe is the only data point corresponding to the peak-effect region 

^ There is some confusion in Ref. [18]. While Fig. 2 provides Hc2 = 50 kOe, the 
value of Hc2 evaluated from Fig. 1 is closer to 45 kOe. Taking into account that 
the resulting A is not very sensitive to some variation of the assumed Hc2 value, we 
have chosen Hc2 = 47.5 kOe for the analysis. 
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Fig. 3. cr (upper panel) and asc = ■\Jcr'^ — cf^g (lower panel) versus H for a CeRu2 
sample studied in [18]. The solid lines represent the theoretical curves calculated 
as explained in the text. The chosen value of Hc2 = 47.5 kOe is indicated in the 
figures by vertical lines. Only data-points shown by closed symbols were used for 
evaluation of ai^g and A. 



(see Fig. 1 in [H]). Because the origin of this effect is not yet established, we 
exclude the corresponding data point from the analysis. 

As may be seen in Fig. 3(a), all data-points for if < 30 kOe can be fairly well 
fitted by the theory, providing A = (167 ± 3) nm and ai,g = (8.4 ± 1) G. The 
magnetic field dependence of age = — a^g is shown in Fig. 3(b). 

Using the asc{H) plot presented in Fig. 3(b), we can calculate A for each of 
the data-points. Such calculations were made for two different values of ifc2 
and they are presented in Fig. 4. As was expected, the absolute value of A 
is practically independent of the chosen value of ifc2- One can also see that, 
contrary to claims of Ref. [E], there is no any noticeable dependence of A on 
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Fig. 4. A as a function of H calculated assuming Hc2 = 45 kOe and Hc2 = 50 kOe. 
H. 

At the same time, the value of a(40kOe) deviates quite significantly from 
the theoretical curve (see Fig. 3)l2J If this deviation is not an experimental 
error, it means that the distribution of the magnetic induction in the case of 
the peak-effect is rather different in comparison with the conventional mixed 
state. However, one should be extremely careful with such conclusions. In the 
case of the peak-effect, the value of a is rather sensitive even to insignificant 
variations of H (see Fig. 13(d) in [18]). In this situation, ~ 1Q~^H change 
of the external field may explain the difference between cr(40kOe) and the 
theoretical curve. 



4-3 Vanadium single crystal 



We discuss experimental data of Ref. [22] in some detail in order to demon- 
strate general problems of interpretation of /iSR experiments in the case of 
low-K superconductors. We also discuss some typical errors that can be found 
in the literature. 

Vanadium is one of the very few pure metals, which displays type-II super- 
conductivity at all temperatures. Superconducting characteristics of vanadium 
have been rather well investigated (see, for instance, [M|35f36ll57] ). Although 
vanadium has a cubic (bcc) structure, Hc2 depends on the orientation of the 
applied magnetic field [35J. According to [35], the value of Hc2 along [111] 
direction is approximately 10% higher than that for [001]. 

^ cT(40kOe) is larger than the corresponding theoretical value. The higher a means 
smaller A. This conclusion is just opposite to that made in the original publication. 



13 



^2 







o 



2- 



0. 



"I r 



(a) 



(b) 



[D 



""I I I I I I 
Vanadium 

□ Ref. 34, if II [110] 
• Ref. 35, if II [110] 

a o Ref. 35,ii||[lll] 

m 

□ 

□ 



H \- 



H \ 



□ Ref. 36, ii II [491] 
E Ref. 37, ii II [491] 



J I L 



J I L 



1 2 3 4 5 

T(K) 

Fig. 5. Hc2{T) for three different orientations of vanadium single crystals 



The dependences Hc2{T) for three different orientations are shown in Fig. 5. 
As may be seen, results of different studies are in excellent agreement. The 
value of Tc may be evaluated as = (5.40 ± 0.05) K [5HI35II36|I57] . While Hc2 
is orientation dependent, its normahzed temperature dependence is practically 
universal [35]. This is illustrated in Fig. 6 where Hc2(T) / (TcdHc2/ dT)T=Tc is 
plotted versus T/T^. The results of different works, presented in such a way, 
nicely collapse onto a single curve. We note that the temperature dependence 
of Hc2 is somewhat different from predictions of Helfand-Werthamer (HW) 
theory [55] . 

Vanadium is a Iow-k material with k{OK) = 1.5 for a pure sample investigated 
in [51]. This circumstance adds some peculiarities to the mixed state and its 
description. First, the condition H <^ Hc2 is not satisfied even in magnetic 
fields down to ifci, i-e., the London approach, in which vortices are considered 
as independent, is inapplicable in the entire range of magnetic fields. In this 
situation, the actual magnetic induction distribution in the sample strongly 
depends on spatial variations of the order parameter, and the accuracy of the 
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Fig. 6. The normalized upper critical field of vanadium single crystals as a function 
of t = T/Tc- The solid line represents the HW theory [38j. 

corresponding calculations plays a crucial role. Second, if /t ~ 1, the condition 
\ ^ (^p = 0.74^(0K) is the Pippard coherence length) is not satisfied 
at low temperatures and electrodynamics become nonlocal, i.e., quantitative 
applicability of the GL theory at T <^ Tc is questionable. Furthermore, the 
results of Refs. [M|35ll36j clearly demonstrate that superconducting properties 
of vanadium at T -C Tc cannot be described by the GL theory and a more 
complex approach is necessary. At the same time, as we argue below, experi- 
mental (t{H) curves are close to theoretical predictions of Brandt [5] and can 
be used for evaluation of the magnetic field penetration depth. 

The values of a were calculated using the P{B) curves presented in Ref. [22]. 
Clearly visible peaks arising from muons stopped outside the sample were 
approximated by Gaussians and subtracted from the data. The resulting values 
of cr are plotted in Fig. 7(a) Our analysis gives A = (49 ±1.5) nm and 
Hc2 = (3.8 it 0.15) kOe. The estimation of A is in very good agreement, 
with A = 50 nm, which may be calculated from Hc2{T) and k,{T) curves 
experimentally measured in 

The value of Hc2 evaluated above is just 10% below of ifc2 = 4.2 kOe provided 
in the original publication [22]. We note that there are no estimations of 
experimental uncertainty for Hc2 in [22]. One can assume that ifc2(0.02K) was 
obtained by extrapolation of higher temperature data and the corresponding 
error margins are considerable. We also note that ifc2(0.02K) = 4.2 kOe is 
well above earlier results (see Fig. 5). Partly this difference may be explained 



^ Because the demagnetizing factor of the sample is close to 1, one can safely assume 
H = B for all considered magnetic fields. 
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Fig. 7. (a) a versus H for a vanadium single crystal studied in [22]. The solid line 
represent a fit with the Brandt theory [5]. The resulting values of A and Hc2 are 
indicated in the figure. 

by the fact that the sample that we are discussing here was substantially less 
pure than those of Refs. [M|35f36ll^ . However, such a significant increase of 
Hc2 seems to be unlikely. Furthermore, as we show below, the value of ifc2, 
evaluated by the analysis of the temperature dependence of muon relaxation 
rates, agrees better with the estimate of Hc2 presented in Fig. 7 than with the 
value given in [22] . 

As was already mentioned, experimental results presented in [3l] allow for 
evaluation both Hc2(T) and A(T) dependences. Using these data, we can also 
obtain the expected value of a for any magnetic field and temperature. Such 
results for if = 1.6 kOe are plotted in Fig. 8 for comparison with the /iSR 
data of Ref. J22|. As may be seen in Fig. 8, the two cr(T) curves are similar. In 
order to emphasize this similarity, we approximate both data sets by the same 
functional dependence (see Fig. 8). Considering the results presented in Fig. 8, 
one can conclude that the sample investigated in [22] has indeed a somewhat 
higher Hc2- The value oiTc{H) may straightforwardly be evaluated from cr(T) 
data as the value of T, at which a vanishes. Such estimate gives Hc2 = 1-6 kOe 
at T = 3.38 k|3] Using this value of ifc2(3.38K) and the normalized Hc2{T) 
curve presented in Fig. 6, we can evaluate Hc2{0) = (3.4 ± 0.25) kOe, which 
is in reasonable agreement with the estimate made from the analysis of cr{H) 
data (see Fig. 7)L^ 

^ This is well below the value of Tc(1.6kOe) = 3.65 K provided in |22j . 
^ We use Tc = 5.4 K, as it follows from earlier measurements (the same value is 
provided in reference data of Goodfellow Ltd.), assuming that Tc = 5.2 K given in 
[22], is a misprint. If, however, we except Tc = 5.2 K, Hc2{0) = (3.65 it 0.2) kOe, 
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Fig. 8. Temperature dependence of cj(1.6kOe) for a vanadium single crystal studied 
in |22j . The temperature dependences of a and A^^ (right y-axis), evaluated from 
experimental data of j34 ]. are shown for comparison. The solid and the dashed lines 
are the guides to the eye. The dotted line (/'(T)) is obtained by scaling of the solid 
line (/(T)), i.e., /'(T) = 1.29/(1.09r). 

Our results presented in this section are rather different from the conclusions 
of Ref. [22] • First and foremost, as it is clearly demonstrated in Fig. 7, the 
magnetic field dependence of a is very close to the result of the GL theory. 
Our value of A for T = 0, which is in excellent agreement with measurements 
of Ref. |34j, is about 1.5 times smaller than the result of p2j for H = 1.6 kOe. 
The temperature dependence of 1/A^ calculated according to Ref. is also 
plotted in Fig. 8. As may be seen, while a vanishes at T = 3.12 K, the value 
of A(3.12K) remains almost the same as at T = 0. In other words, cr vanishes 
at Tc{H) not because of the divergence of A but because the coefficient F in 
Eq. (7) vanishes at this temperature. 

It seems important to emphasize that Fig. 10 of Ref. [22] is based on an obvious 
misunderstanding^ There exists no theory that predicts divergence of A at 
Tc{H). The reference on theoretical calculations of Miihlschlegel [39], given 
in [22] , is misleading. Indeed, the thermodynamical consideration of Ref. [SH] 
is based on the fact that the difference between free energies of the normal 
and the superconducting states per unit of volume can be written as H^/Stc. 
The same difference can also be written as ricpA where ricp is be the density 
of Cooper pairs and A is the equilibrium (zero-field) superconducting energy 
gap. Using this, one obtains UcpiT) and A(T). Nothing in this consideration 
can be used to justify Fig. 10 of Ref. [22]. We also note that the temperature 
variation of Hc2 should be taken into account if the temperature dependence of 
A is evaluated from measurements in fixed magnetic fields. It can be neglected 



the value practically coinciding with the result of Fig. 7. 
^ A similar plot one can also find in [30j 
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only if the condition H ^ Hc2 is satisfied at all temperatures. 



In fact, good agreement with the theory, demonstrated in Figs. 7 and 8, is 
rather surprising. As was already mentioned, it is well established that vana- 
dium does not obey the GL theory at T <^ Tc [Mll35ll36] .The most probable 
is that the distribution of the magnetic induction in the sample {P{B)) is 
different from theoretical predictions, while a, as a more integral character- 
istic of this distribution, remains practically the same. This assumption can 
also explain the difference between our results and those of Ref . [22] . Indeed, 
the analysis of P{B) functions, carried out in [22j, resulted in an unphysical 
magnetic field dependence of A, which clearly demonstrates the inapplicability 
of the GL theory to this analysis. 

It is important to underline that, although the distribution of the magnetic 
induction in the sample cannot be described by the GL theory in Iow-k type-II 
superconductors at T <^ Tc, the cr{H) curves can still be used for evaluation of 
the magnetic field penetration depth, as it is proven by a very close agreement 
between our value of A (OK) = 49 ib 1.5 nm and 50 nm calculated from results 
of Ref. [31]. We also note that at temperatures closer to Tc the GL theory 
should be applicable and both analyses should result in the same A values. 



44 YNi2B2ClT]^ 



In order to demonstrate that in some cases the GL theory cannot describe 
/iSR data, we consider a study of a borocarbide superconductor YNi2B2C |T^ . 
Rare-earth nickel borocarbide superconductors attracted a lot of attention dur- 
ing the past decade. Already the very first studies of YNi2B2C demonstrated a 
pronounced positive curvature of the Hc2{T) curve, indicating unconventional 
superconductivity jlHllI]. Similar conclusions were made from specific heat 
data pOlfiT] . Although YNi2B2C has been extensively studied, the nature of 
this unconventionality is still under discussion. While Refs. |42|43|44|45|46II47] 
provide evidences of point nodes in the superconducting gap function, other 
works point out on multiband superconductivity [48,49,50.51] . The distinction 
between these two possibilities is sometimes difficult to make. For instance, as 
was recently demonstrated, specific heat data may be fitted equally well by 
nodal and two-gap models [52]- 

Experimental results of Ref. [14] are plotted in Fig. 9 as a versus H. The 
value of ifc2(3K) = 70 kOe for this particular crystal is given in As may 
clearly be seen in Fig. 9, cr{H) data-points cannot be fitted with the theory if 
the entire range of magnetic fields is considered. Because there are sufficient 
experimental evidences that YNi2B2C is an unconventional superconductor 
[,42,4^pill45f4tif47ll48f49f50ll51f52] . disagreement between the GL theory and 
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Fig. 9. (T as a function of H for a YNi2B2C sample studied in Ref. ^3]. The solid 
lines are the theoretical a{H) dependencies calculated for Hc2 = 70 kOe and for 
two different values of A. 

experimental data is expected. We also note that in YNi2B2C a transition 
from a triangular to a square vortex lattice was observed [53)151] . However, 
because this transition occurs in lower r nag netic fields, it cannot have any 
influence on /xSR data presented in Fig. qLII 

While the totality of data cannot be fitted with the theory, both high-field and 
low-fild results may amazingly well be approximated with two different theo- 
retical curves, corresponding to two different A values (Fig. 9). Unfortunately, 
insufficient number of data-points does not allow to make unambiguous con- 
clusions on this matter, however, if this behavior will be confirmed by a more 
detailed study, it may be considered as a rather interesting result, indicating 
two gap superconductivity. 

In low magnetic fields H <^ ifc2, niost of muons stopped outside vortex cores, 
i.e., the magnetic induction distribution in vortex core regions is not very 
important for /iSR data. In this case, the difference between conventional and 
two gap superconductors should not be significant and the resulting a{H) 
curves can be close in these two cases. 

In higher magnetic fields, as it was established in studies of MgB2, supercon- 
ductivity in one of two bands is completely suppressed and the superconduc- 
tor behaves itself as a one gape superconductor but with a smaller number of 
Cooper pairs [55|56II57] . This can explain the fact that the two data points for 




^ As was demonstrated in [10], although the magnetic induction distributions for 
square and triangular lattices are quite different, <j{H) remains practically the same 
in both cases. 
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H > 30 kOe follow a standard theoretical curve with a higher value of A (see 
Fig. 9). 

The quantity 1/ is proportional to the density of Cooper pairs. If two gap 
superconductivity is assumed, the values of the magnetic field penetration 
lengths, evaluated from low- and high-field data, allows evaluation of relative 
weights of two superconducting bands. Such estimate gives 54% and 46% for 
stronger and weaker gaps, respectively. These values are noticeably different 
from the result 71% and 29% obtained in Ref. [52]. At present, however, it is 
too early to discuss such differences. Two data points in the high-field range 
part of the curve (see Fig. 9) are clearly insufficient in order to make any 
definite conclusion about superconductivity in YNi2B2C. 



5 Conclusion. 

In this work, we applied numerical calculations of Brandt [5] for the analysis 
of /iSR experiments carried out in the mixed state of several superconducting 
compounds. It turned out that this approach may serve as a very powerful tool 
for the interpretation of /xSR experiments. If the magnetic field dependences 
of muon depolarization rates are available, not only A but also Hc2 can reliably 
be evaluated. We show that in the most of considered cases the magnetic field 
dependences of a may very well be described by a single and temperature 
independent A. 

In contrast to approximate analytical models, Ref. [5] provides precise numer- 
ical solutions of 2-dimensional GL equations for different values of k (0.85 < 
K < 200) and for magnetic fields ranging from Hd to Hc2- Using these so- 
lutions, different characteristics of the mixed state, including the a{B/Bc2) 
dependences for various k values, were calculated. As well as we are aware, 
these calculations provide the best available description of the magnetic induc- 
tion distribution in the mixed state of conventional type-II superconductors. 
We also note that numerical calculations of cr{H) are available since 1997 (see 
Fig. 3 in [1]). In spite of this, for some mysterious reasons, these theoretical 
calculations have practically never been used for the analysis of /iSR data. 

We also argued that the magnetic field dependence of A can never be obtained 
from analyses of experimental data collected in the mixed state. Indeed, be- 
cause the local value of A is inversely proportional to the absolute value of the 
superconducting order parameter, one cannot introduce any single value of A 
in the mixed state. 

Calculations of Brandt [5] represent the conventional GL theory and their 
validity for the description of unconventional superconductors is question- 
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able. In fact, there are no reasons to believe that the conventional GL theory 
can quantitatively describe either two-gap superconductors or superconduc- 
tors with nodes of the order parameter and one should expect disagreement 
between Brandt's theory and experimental results in the case of unconven- 
tional superconductors, as it is demonstrated in Fig. 9. 

We demonstrated that in conventional superconductors, the results of /xSR 
experiments may be used for the evaluation of both A and Hc2- If applicability 
of the conventional GL theory is questionable, the knowledge of Hc2 is of 
primary importance. Disagreement between the values of resulting from 
/iSR data and independent measurements may be considered as convincing 
evidence that this particular superconductor is unconventional. 

We are grateful to R. Khasanov for numerous and fruitful discussions. 
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